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Abstract. With a peak luminosity of ~ 10 47 erg s , the December 27th 2004 giant flare from SGR1806-20 would have been 
visible by BATSE (the Burst and Transient Source Experiment) out to ~ 50 MpcQLQ]. It is thus plausible that some fraction 
of the short duration Gamma-Ray Bursts (sGRBs) in the BATSE catalogue were due to extragalactic magnetar giant flares. 
According to the most widely accepted current models, the remaining BATSE sGRBs were most likely produced by compact 
object (neutron star-neutron star or neutron star-black hole) mergers with intrinsically higher luminosities Q]. Previously, by 
examining correlations on the sky between BATSE sGRBs and galaxies within 155 Mpc, we placed limits on the proportion of 
nearby sGRBs (3]. Here, we examine the redshift distribution of sGRBs produced by assuming both one and two populations 
of progenitor with separate Luminosity Functions (LFs). Using the local Galactic SGR giant flare rate and theoretical NS-NS 
merger rates evolved according to well-known Star Formation Rate parameterisations, we constrain the predicted distributions 
by BATSE sGRB overall number counts. We show that only a dual population consisting of both SGR giant flares and NS-NS 
mergers can reproduce the likely local distribution of sGRBs as well as the overall number counts. In addition, the best fit LF 
parameters of both sub-populations are in good agreement with observed luminosities. 
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INTRODUCTION 

The leading candidate progenitor model for short dura- 
tion Gamma-Ray Bursts (sGRBs) is the merger of two 
compact objects, neutron star-neutron star (NS-NS) or 
neutron star-black hole [1]]. However the initial spike 
in a giant flare from a Soft Gamma Repeater (SGR) 
in a relatively nearby galaxy would also appear as a 
sGRB. For example, the December 27th 2004 event from 
SGR 1806-20 would have been visible by BATSE out to 
fa 50 Mpc fH H. Thus it is plausible that some frac- 
tion of sGRBs are extragalactic SGR giant flares. Pre- 
viously, Tanvir et al. [3] demonstrated that between 10 
and 25 per cent of BATSE sGRBs were correlated on 
the sky with galaxies within « 113 Mpc. The Luminos- 
ity Function (LF) of sGRBs has been investigated previ- 
ously assuming a single population [4] in order to deter- 
mine the intrinsic sGRB rate and most likely LF param- 
eters. Here we assume intrinsic rates given by both the 
observed Galactic SGR flare rates and modelled NS-NS 
merger rates in order to investigate single and dual pop- 
ulation LF parameters. Obviously there are significant 
uncertainties in these rates: the Galactic giant flare rate 
in particular is estimated from only 3 observed events. 
Regardless of these uncertainties and the exact form of 
LF for a second population, we find a single progenitor 
population constrained by overall number counts cannot 
produce sufficient local events. We present preliminary 
results that show only a dual population can reproduce 



the likely local sGRB distribution as well as the overall 
number counts. 



METHODS 

The number of sGRBs, N, observed above threshold p in 
time T and solid angle Q. is given by equation Q] where 
<3>(L) is the sGRB LF, Rgrb{z) is the sGRB event rate at 
redshift z, dV (z) / dz is the comoving volume element at 
z and Zmax f° r a burst of luminosity L is determined by 
the detector flux threshold and the luminosity distance 
of the event. We assume a spectral photon index of 1 . 1 in 
the BATSE energy range to convert photon flux to energy 
flux. 
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(1) 

The sGRB rate, Rgrb(z) i s gi ven by equation [2] where 
Pgrb is the rate of sGRBs per progenitor, pp(0) is the 
intrinsic progenitor formation rate and F (z) describes the 
progenitor production rate as a function of z- 



Rgrb = Pgrb x Pp(0) x F(z) 



(2) 



For NS-NS mergers, Pgrb = 1 and is the delayed 
rate given by the convolution of the Star Formation Rate 
at z (SFR(z)) with a distribution of delay times to repre- 



sent merger formation. For SGR flares, Pgrb = 3x 10 
(the observed local rate of giant flares per Galactic SGR) 
and F(z) follows both SFR(z) for magnetar production 
from supernovae and the delayed SFR(z) to allow for 
production by White Dwarf binary mergers |5|]. We take 
SFR(z) as given by the SF2 model of Porciani and Madau 
|@], with the local SFR from 0]. P/»(o) SGR rates per 
Milky Way equivalent galaxy are from Levan et al. |5|], 
and we follow Kalogera et al. [8] for a NS-NS merger 
rate of 10~ 5 yr _1 . 

LFs for SGR giant flares and NS-NS mergers are 
not well constrained. A log-normal LF approximates the 
theoretical NS-NS merger distribution [9], and is also 
plausible for SGR giant flares. Here we assume sGRBs 
(of differing luminosities) to arise from both SGRs and 
NS-NS mergers. We also investigate a single merger 
population with Schechter function and log-normal LFs. 
Merger time distributions are similarly unknown, and we 
assume P(log(T))dlog(r) = constant between 10 and 
10 10 years, approximating theoretical distributions 1 10]. 

The Cmax/Cmin table from the BATSE catalogue ltllll 
provides peak count rate for bursts in units of threshold 
rate. The BATSE threshold was varied historically, and to 
analyse a consistent set of bursts we restricted the table 
to only those sGRBs recorded when the threshold was 
set to 5.5<T above background in at least 2 detectors in 
the 50 — 300 keV range. The all sky equivalent period 
this represents is ~ 1.3 years. By varying the parameters 
of the LFs (median and sigma for log-normal, slope 
and for Schechter) and fitting the predicted N(> 
p) to the Cmax/Cmin distribution, we found the best fit 
LF parameters by % 2 minimization. Figure Q] shows an 
example of such a fit. Having obtained the best fit LF 
parameters, the redshift distribution of predicted sGRBs 
can be calculated, as shown in Figures |5]|7] 
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FIGURE 1. Example of best fit plot of number counts pre- 
dicted from Equation Q] using dual population lognormal LFs 
versus observed C max /C mm - Fits are performed for differential 
number counts, but shown cumulative in the figure. Dashed 
line is observed C max /C m i n , solid line predicted number counts. 
Reduced ^2 = 1-58. 
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FIGURE 2. Dual population best fit log-normal Luminos- 
ity Functions from Figure Q] Dotted line is LF for bursts 
due to SGR giant flares, median log(E p ) = 44.5, (7/ og = 0.95. 
Dashed line is LF for bursts due to NS-NS mergers log(E p ) = 
48.0, G/og = 0.95. Other details as Figure[TJ 



RESULTS AND CONCLUSIONS 

Although a single NS-NS merger population luminosity 
function can produce reasonable fits to the overall num- 
ber counts, only samples with two populations of dif- 
fering luminosity distributions reproduce the likely local 
population of sGRBs (Figures [3] and [5]l - even Schechter 
type single luminosity functions, which are dominated 
by low-luminosity events, fail to reproduce the nearby 
population (Figure [7J- Furthermore, the overall redshift 
distribution of the dual population samples are consis- 
tent with the handful of redshifts known for sGRBs lfL?!] . 
Despite large uncertainties in assumed intrinsic rates, the 
best fit LF parameters also agree well with the few peak 
luminosities currently known for both SGR giant flares 
and cosmological sGRBs. 
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FIGURE 3. Cumulative distribution of bursts within 500Mpc 
using best fit dual population log-normal LFs results from 
Figure [2] Also shown (open squares with lo" errors) are the 
rates of local sGRBs as measured by correlation analyses with 
local galaxies from Tanvir et al. Q]. Other details as Figure[TJ 




FIGURE 4. Cumulative distribution of bursts out to z = 3 
using best fit dual population log-normal LFs results from 
Figure[2| Also shown (solid line) is the cumulative distribution 
of the handful of sGRBs with known redshifts |12]. Other 
details as Figure[T] 
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FIGURE 7. Cumulative distribution of bursts within 500Mpc 
using best fit single Schechter function (L* = 50.25, slope = 
—0.9, reduced %2 — 2.53.) for the NS-NS merger population. 
Other details as Figure[3] 
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FIGURE 5. Cumulative distribution of bursts within 500Mpc 
using best fit dual population with log-normal LF for the SGR 
flares (median log(E p ) = 45.0, 0/ og = 0.75), and a Schechter 
function (L, = 50.0, slope = -0.9) for the NS-NS merger 
population. Reduced X2 = 2.36, other details as Figure[TJ 
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FIGURE 6. Cumulative distribution of bursts out to z = 
3 using best fit dual population (lognormal SGR, Schechter 
mergers) as Figure[5] Other details as Figure|4] 
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